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1 Introduction 



The notion of the reflection equation associated with solutions of the Yang-Baxter equation PJ |2] , 
goes back to the key works of Cherednik |3] and Sklyanin jl]. The subject has recently attracted a 
great deal of activity as was summarised in [Hj (and references therein) . More specifically, starting 
from a quantum R- matrix -R(A) depending on the spectral parameter A and satisfying the (super) 
Yang-Baxter equation [TJ El 

Ri 2 {Xi - A a ) i2i 3 (Ai) i? 23 (A 2 ) = # 23 (A 2 ) i?i 3 (Ai) i2ia(A a - A 2 ) , (1.1) 

one derives the reflection equation for an object K(X) as 

i?i 2 (Ai - A 2 ) ^i(Ai) R 12 (\i + A 2 ) K 2 (X 2 ) = ^ 2 (A 2 ) R 12 (Xx + A 2 ) ^(Ai) fl 12 (Ai - A 2 ) . (1.2) 

We have proposed in [5] a classification of c-number solutions if (A) of the reflection equation (jl.2j) 
for rational (super) Yangian R- matrices |H] associated to the infinite series so(m), sp(2n) and 
osp(m\2n) . This classification entailed K matrices with purely diagonal, anti-diagonal and mixed 
(diagonal, anti-diagonal) non-zero entries. The explicit values of the K matrices were then used 
within the analytical Bethe Ansatz formulation [HI E3 El El Ej for the derivation of the spectrum 
and the bulk and boundary ^-matrices for the so(m), sp(2n) open spin chains. 
There exists a substantial body of work on gl(m\n) super spin chains. Interest in these systems 
stemmed from the existence of physically relevant particular cases such as supersymmetric t-J 
and extended Hubbard models. They have been the object of many studies. Supersymmetric 
t-J models were considered, e.g. in (thermodynamical aspects), ^1] (diagonal boundary K 
matrices) and ^1 (boundary S'-matrix). Extended Hubbard models were considered in JE] (closed 
chain) and in [13 CHI CHI (open chains with integrable boundary conditions), whilst spin ladder 
systems associated to some sl(m\n) superalgebras were obtained in (201 [2T]. General results for 
continuum limit of the gl(m\n) super spin chains were derived in |22j . 

A natural alternative to these models with gl(m\ri) underlying superalgebras is provided by super 
spin chains with underlying osp(m\2n) superalgebras. A connection to intersecting loop models 
and hence polymer field theories was pointed out in [2S], where the analytical Bethe Ansatz 
equations were written for the closed spin chain. Algebraic methods were used for some specific 
cases in |24|l25] including nested Bethe Ansatz in j2Hj. Field theoretical limits were also considered 
in the literature: the exact bulk osp(2\2) S'-matrix was conjectured in [2B] in the framework of 
disordered systems. Investigation of the thermodynamics of osp(l\2n) closed spin chains was 
undertaken in [2Z| using the thermodynamical Bethe Ansatz formalism. An algebraic construction 
using Birman-Wenzl-Murakami algebra then yielded conjectural S-matrices for field theoretical 
osp(m\2n) models, and allowed a subsequent thermodynamical Bethe Ansatz analysis of their 
thermodynamical properties [20] . However a systematic thermodynamic treatment of these models 
with more general boundaries is still missing. 

Our purpose is to make an exhaustive study of the more complicated case of open spin chains 
with osp(m\2n) underlying superalgebra and any integrable (diagonal at a first step) boundary 
conditions. The strategy is to establish (insofar as the methods are available) Bethe Ansatz 
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equations for ground state and excited states (note that there is no obvious relation between closed 
spin chains Bethe Ansatz equations and open spin chain Bethe Ansatz equations, particularly 
when non trivial boundary conditions are involved); solve them within the non trivial string 
hypothesis (discussed in the closed case in [THJ 123); and use the results to obtain the ^-matrix 
and thermodynamical quantities, with explicit evaluation of the effect of boundary conditions. 

This paper is our first step in this direction: using the analytical Bethe Ansatz method, we 
derive the Bethe Ansatz equations for all orthosymplectic superalgebras, and all diagonal K 
matrices. Restricting ourselves then to osp(l\2n), we solve these equations in the thermodynamic 
limit, we derive the ground state and low- lying excitations, and compute explicitly the bulk and 
boundary S-matrices. Further generalisations will be left for future investigations. 

2 Bethe Ansatz equations for the osp(M|2n) open spin 
chain 

2.1 Conventions and notations 

The Bethe Ansatz equations will be derived here for the osp(M |2n) iV-site open spin chain with 
diagonal reflection conditions by means of the analytical Bethe Ansatz method (see e.g. [3 EH 02] )■ 
As customary to construct the open chain transfer matrix we introduce the i?-matrix which is a 
solution of the super Yang-Baxter equation. We focus on the osp(M\2n) invariant .R-matrix given 
by i 

R(\) = \(\ + i K )l + i(\ + iK)P-i\Q , 2k = 6 (M - 2n - 2) (2.1) 
where P is the (super)permutation operator (i.e. X 2 i = PXi 2 P) 

M +2n 

P= J2 (-l) b '%®^ (2-2) 

and 

M+2n 

Q = (-l^OiOjEjt ® Eji = P* 1 . (2.3) 

For each index i, we have introduced a conjugate index 

t = M + 2n + l-i. (2.4) 

We also introduce a sign 6^ and a Z 2 -grading [i] whose definition, due to the conventions we adopt 
(see below), depend whether we consider the superalgebra osp{2\2n) or any other osp(M\2n) 
superalgebra: 

For osp(M\2n) superalgebras, M^2: 

,-1 for l<i<M + n 

■1 for M + n+l<i<M + 2n 
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(_1)M=+1 for l<i<n and M + 1 < i < M + 2n (2.6) 
(-1)W = -1 for ri + 1 <i <n + M (2.7) 

We will associate to this choice the sign # = — 1- 
For osp(2\2n) superalgebras: 

f+1 for l<z<n+l and % = 2n + 2 . . 

Oi = < ~ ~ 2.8 

[-1 for n + 2<i<2n + l 

(-l)M = +1 for z = 1 and z = 2n + 2 (2.9) 
(-l)M = -1 for 2 < i < 2n + 1 (2.10) 

The sign corresponding to this choice will be 8$ = +1. 

The transposition 4 used in (12. 3|) and below is defined, for A = . A y E^, by 

A* = ^(-1)M[?1+[J]^ A iJ £ jT = ^ (A*) ij Eij (2.11) 

The i?-matrix (j2.1|) satisfies crossing and unitarity, namely 

i?i 2 (A)i?i 2 (-A) = (A 2 + K 2 )(A 2 + 1)1, i? 12 (A) = i?i 1 2 (-A-m). (2.12) 
We finally define the super trace operation according to the Z 2 -grading we have introduced: 

M+2n M+2n 
j=l i,j=l 

2.2 Transfer matrix and pseudo vacuum 

The open chain transfer matrix is given by 0] 

t(A) = Tr K+(A) T (A) Kq(\) f (A) , (2.14) 

where Tr denotes here the super trace ()2.13|) over the auxiliary space, 

T (A) = R QN (X)Ro, N ^(X) ■ • ■ R 02 (X)Roi(X) , T (A) = R W (X)R 20 (X) ■ ■ ■ R N -i, (X)R N0 (X) , (2.15) 

Kq(X) is any solution of the super boundary Yang-Baxter equation 

R 12 (Xt - A 2 ) K^Xx) R 12 (Xt + A 2 ) K 2 (X 2 ) = K 2 (X 2 ) R 12 {X X + A 2 ) K 1 (X 1 ) R 12 (Xt - A 2 ) (2.16) 

and Kq (A) is a solution of a closely related reflection equation defined to be: 

-Ri 2 (A 2 -Ai) Kl^Xt) R l2 {-Xx-X 2 -2m) K%{X 2 ) = K%{X 2 ) R 12 {-X x -X 2 -2iK) K£(\ x ) R 12 (X 2 -X 1 ) . 

(2.17) 



It is clear that any solution K (A) of (j2.16|) . e.g. given in [S], gives rise to a solution K + {\) of 
dHZJ), defined by K+(X) = K~(-X - inf. 

To determine the eigenvalues of the transfer matrix and the corresponding Bethe Ansatz equa- 
tions, we use the analytical Bethe Ansatz method [3 EH E] We follow the same procedure as in 
by imposing certain constraints on the eigenvalues, deduced from the crossing symmetry of 
the model, the symmetry of the transfer matrix, the analyticity of the eigenvalues, and the fusion 
procedure for open spin chains. These constraints allow to determine the eigenvalues by solving a 
set of coupled non-linear consistency equations or Bethe Ansatz equations. 

We first describe the case with trivial boundaries, K~(X) = K + (\) = 1. 

We recall that the fusion procedure for the open spin chain ^21 I2H| yields the fused transfer matrix 
t(A) = C(2A + 2ik) t(X) t(X + zk) - C(A + m) 2Af g(2A + w)q(-2X - 3m) , (2.18) 
where we define 

C(A) = (A + m)(A + i)(X- ik)(X - i) , q{X) = 6 (X - i6 )(X - we). (2.19) 

Note that the value of g(A) is related to the specific choice of the position of the orthogonal and 
the symplectic part in the .R-matrix. We choose for the general case (apart from the osp(2\2n) 
case) the symplectic part to be "outside" and the orthogonal part to be "inside" . This formulation 
corresponds to a specific Dynkin diagram: the so-called distinguished one (see fig. HJ). For the 
case of osp(2\2n) the distinguished Dynkin diagram has a special form (see fig. Q): it corresponds 
to the orthogonal part being "outside" and the symplectic part being inside. These considerations 
justify the conventions we have adopted in ()2.5|) - (j2.6|) and ()2.8|) - (j2.9|) . 

From the crossing symmetry of the i?-matrix (J2.12j) it follows that: t(A) = t(— X— ik). The transfer 
matrix with K~ = K + = 1 is obviously osp(M\2n) invariant, since the corresponding i?-matrix 
(|2.1j) is osp(M\2n) invariant, namely 



0, (2.20) 



where U is any generator of the osp(M\2n) algebra. Finally, from the assumption of analyt- 
icity of the eigenvalues, we require that no singularity appears in the Bethe eigenvalues. The 
aforementioned set of constraints uniquely fix the eigenvalues. 

We now choose an appropriate pseudo-vacuum, which is an exact eigenstate of the transfer 
matrix; it is the state with all "spins" up, i.e. 



N 

= (R) wnere 1+) 



i=l 



/ 1 \ 





V / 



G C M+2n . (2.21) 



Our choice of #0 ensures that this state is always bosonic, whichever orthosymplectic superalgebra 
we consider. 
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O" o ^ 

osp{2m + l\2n), m > 1 

o— o — — 0-04=0 



osp(l\2n) 
ai a n -i dr. 

O— - 0=4* 



osp{2\2n) 

d\ 0,2 a n Ctn+l 



Figure 1: Distinguished Dynkin diagrams of the osp(M\2n) superalgebras. 

After some lengthy computation, we determine explicitly the action of the transfer matrix as 
t(X)\u+) = A°(A)|c<o)) where A°(A) is given by the following expression 

2n+M-2 

A°(A)=a(A) 2 %(A) + 6(A) 2Ar £ (-l) [l+1] 9l (X) + c(A) 2 % n+M -i(A) (2.22) 

1=1 

with 

a(\) = (\ + i)(\ + iK), 6(A) = A(A + i/c), c(A) = A(A + in - i) (2.23) 
The expressions of the functions gi(X) depend on the case we consider. 
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For the generic osp(M\2n) case, M^2, they are given by (with M = 2m or M = 2m + 1): 

m A(A + y - |)(A + in) 

QiW = : r tttt; ; t = 0, . ..,n— 1, 

(A + f)(A + f)(A-* (m) ^ 



/u A(A + y - |)(A + in) 

qAX) = : , / = n, . . . , n + m — 1 

i; (A + f)(A + m-f)(A + m-^)' 

3n+ m (A) = A(A , + m) ifM = 2m + l 

y + v ; (A + i^)(A + z^f±±) 

^(A) = (? 2 n + M-n(-A-m), / = 0,l,...,M + 2n (2.24) 

In this case, we have k = n + 1 — 4^. We also set = n + m. 



For the case of osp(2\2n), due to the different conventions, one has: 



(A + f + i)(A + i«) 



(A + f)(A + |) 

A(A + f + i)(A + »«) 

Pa A) = -f Tjj— , f = l,...,n, 

(A+f)(A + i-f)(A + i-^)' 

9iW = «an+i-i(-A-iK), Z = 0,...,2n+1 (2.25) 



We remind that in this latter case, k = — n. 



2.3 Dressing functions 



From the exact expression for the pseudo-vacuum eigenvalue, we introduce the following assump- 
tion for the structure of the general eigenvalues: 

2n+M-2 

A(A) = a(X) 2N g (X)A (X) + b(X) 2N ]T (-1)V + \(X)MX) + c(A) 2 % n+M -i(A)A 2n+M _ 1 (A) 

i=i 

(2.26) 

where the so-called "dressing functions" Aj(A) need now to be determined. 
We immediately get from the crossing symmetry of the transfer matrix: 

Ai(X) = A 2n +M-i-i(—X — in) l = 0,...,M + 2n-l. (2.27) 

Moreover, we obtain from the fusion relation ()2.18|) the following identity, by a comparison of the 
forms (|2.26|) for the initial and fused auxiliary spaces: 

A (X + iK)A 2n+M -i{X) = 1 . (2.28) 

Gathering the above two equations (|2.27|) . ()2.28|) we conclude 

A (A)A)(-A) = 1. (2.29) 
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Additional constraints are then imposed on the "dressing functions" from analyticity properties. 
Studying carefully the common poles of successive gi's, we deduce from the form of the gi functions 
(I2.24|) that g\ and have common poles at A = — ^ or A = —in + |, therefore from analyticity 
requirements 

M~) = A^(~), i = l,...,n-l, 

il il 
A i(-in + -) = A l _ 1 {-in + -), I = n, ...,n + m-l (2.30) 

There is an extra constraint when M = 2m + 1, namely 

A n+m (-in + y) = A n+m _i(-m + y) . (2.31) 

Having deduced the necessary constraints for the "dressing functions", we determine them ex- 
plicitly. The "dressing functions" A\ are basically characterised by a set of parameters A^ 
with j = 1,...,M^, where the integer numbers are related to the diagonal generators 
of osp(M\2n). These generators are defined as: 

N 

S {l) = J2 s i\ s® = (e u -eu)/2, {e kl ) l3 = 5 ik 5 jt . (2.32) 

The precise identification of follows from the symmetry of the transfer matrix (see also jU]): 
S (i) = M {l ~ 1} -M (l \ Z = l,...,n-l,n+l,...,n + m-2 (2.33) 

S {n) = M (n-1) _ 2M (n) ^ 
s (n+m-l) _ M (n+m-2) _ ^(n+m) _ ^(n+m-l) _ M (n+m) ^ jf ]Vf = 2m + 1 (2.35) 

5(n+m _i) = M {n+m-2) _ M {+) _ M (-) ] 5 (n+m) = M (+) -M^, ifM = 2m (2.36) 
and AfW = f . 

A. osp(2m+ l|2n) 

The dressing functions take the form: 

M<1) A + A« - i A - A?> 



A (A) = n 



J 2 ' v "j 2 



; A + \f + | A- Af } + 



j=l " 1 1 2 " "J" 1 2 



M A + Af + f A-A<" + f 



X 



n A + A j + 2 ~ 2 A ~ A j + 2 ~ 2 , = 1 

j=l A + A j + 2 + 2 A ~ A i +2 + 2 
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Mm \ + \f>+in-%-i A-Af +m-f 



MX) 3 A + A (0 +m _ f A _ A C0 +in _| 

^A + Af^+m-f + l A-A^+m-f + | 

x a- n , Z = n,...,n + m- l 

M A + Af +1) + m-f-i A-A^+m-f-l 

^ A + Af+m-f + i A-Af +m-f + t 
An+m(A) " A + Af + m-f A-Af +m-f 

x A + Af+m-f -f A-Af+m-f -j 
A + Af +m-f + | A-Af+m-f + |' 

and Aj(A) = A 2n+ 2m-i( _ A - m) for / > n + m, « = n — m + \. 



B. osp(2m|2n) with m > 1 

The dressing functions are the same as in the previous case for I — 0, . . . , n + m — 3, but 



X 



X 



X 



i= \ A + Af _2) + in - f + i A - Af " 2) + m - f + i 



m(+) 

n 

M(") 

n 

m(+) 

Wm-l(A) = J| 



A + A™ 


+ in — 


ik 
2 


+ 


3i 
2 


A 


- Aj +) + in 


ik 
2 


1 3« 
2 


A + A« 


+ in — 


ik 
2 


+ 


i 
2 


A 


- Aj +) + %n 


ik 
2 


~ 2 


A + A<"> 


+ in — 


2 


+ 


3i 
2 


A 


- \ { ~ ] + in 


life 
2 


~ 2 


a + a<-> 


+ in — 


ife 
2 


+ 


i 
2 


A 


- aJ _) + in 


ik 
2 


~ 2 


A + Aj +) 


+ in — 


ifc 
2 


+ 


3i 
2 


A 


- Aj +) + in 


ik 
2 


1 3i 
2 


a + a« 


+ in — 


ik 
2 


+ 


2 


A 


- Aj +) + in 


ik 
2 


~ 2 


A + A<"> 


+ in — 


ifc 
2 




i 
2 


A - 


- aJ _) + in - 


2 


j 
2 


A + A" 


+ in — 


ifc 
2 


+ 


i 
2 


A- 


- A^ } + in - 


ifc 
2 


4- - 

T 2 



M(") 

n /v I x> j r ""' 2 2 x> x> j r "'" 2 2 / oo\ 

x , TF) ; z — ^ — x(-) ; z — (2 - 38) 



and Ai (A) = A 2n +2m-z-i(— A — m) for / > n + m — 1, k = n + 1 — m. 



C. osp(2|2n) 

As already mentioned, this case must be treated separately, because of the different position of 
the orthogonal and symplectic parts in the i?-matrix. Here the orthogonal part of the i?-matrix 
is considered to be "outside" and the symplectic part "inside" as opposed to the previous cases. 
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The corresponding dressing functions have the form 



M(i) , . x (i) i x x(i) 



^o(A) = n 



A + A • — 4 A — A ■ — 



2 2 3 2 



J=l A + A j + 1 A - A j + f 

mW , ,(z) i/ , ,(0 jj 

MX) = U J 1 1 



X 



\ A + Af + i - f A - Af + % - § 

M(i+1) \ J_ _L 3i il \ , 3i ii 

n A I A i I I ~ 2 A ~ A j + T ~ 2 7 -, , 

j=l A + A j + 2 ~ 2 A ~~ A j + 2 _ 2 



M(") \ i \( n ) in \ \W in 

n/A /A ■ — TT — — TT 



X 



A + Af +1) +2i-j A - Af +1) +2i-j 

11 x , x(n+l) in x x(n+l) i„ ' ^- 39 ) 

i=i A + A i 2" A ~ A j 2 



and Ai (A) = A 2n+ i-i(— A — m) for Z > n, k — —n. 



2.4 Bethe Ansatz equations 

We define the function 



A 4- — 

e,( A ) = 3— £ • (2-40) 

2~ 

From the analyticity requirements one obtains the Bethe Ansatz equations which read as: 
A. osp(2m+ l|2n) 

mw m( 2 ) 

e^'r = n ^ - x ?) ^ + x f) n «-i(\ w - >f > «-i(\ w + *s*> • 

1 = II e 2 (Af-Af)e 2 (Af + Af) J] II e-^Af - A?^) e_ l( Af + A^) 

j=l,j^i T=±l jr' = l 

/ = 2, . . . , n + m — 1, l^n 
1 = Y[ ei (Af ) -Af + Vi( A S n) + A ? +1) ) II e-i(^ n) -Af-V-i(\- n) + A; n_1) ) 

1 = J] ei (A; n+m) - Af +m) ) ei (AS n+m) + A5 n+m) ) 

j^(n+m-l) 

x II e -i( A !" +m) - A5 n+m - 1} ) e_!(A; n+m) + Aj n+m_1) ) (2.41) 
i=i 
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In particular for M = 1 the Bethe Ansatz equations become 

M W m( 2 ) 



= n ^ - ^ ^ + ^ n e -^ ] - A f ) e -^ ] + ^ 



lf Ji x (2). 

Af« MOW 

1 = II ^(^-Af)e 2 (Af + Af) J] II M^-^MA^+A^) 

j=l,j^i t=±1 j=l 

Z = 2,...,n-1, 

1 = J] e.^A^-Af^e-^A^ + AfVsCA^-Af^e^A^ + Af) 

x II ^(Af - Af- 1} ) e^CAf ) + Af- 1} ) (2.42) 
i=i 

B. osp(2m|2n) with m > 1 

The first n + m — 3 equations are the same as in the previous case for M = 2m + 1, see eq. (|2.41|) . 
but the last three equations are modified, and they become identical to the last three equations 
of the so(2n + 2m) open spin chain, namely, 

1 = H e 2 (X ( t +m - 2) - \ { ; +m - 2) ) e 2 (Af +m - 2) + Af +m - 2) ) 

x II e -i(AS n+m ~ 2) - Af +m " 3) ) e_!(Af +m - 2) + Af +m ' 3) ) 
i=i 

x n J] e^Af^ - Af ) e^Af +m " 2) + A^ } ) 

1 = II e 2 (Af } - Af } ) e 2 (Af } + Af ) J] e_ x (Af } - A< B +m - 2) ) e^A^ + \f +w ~ 2) ) 

(2.43) 

C. osp(2|2n) 

M( 2 ) 

i=i 

1= II e2 (Af-Af)e 2 (Af +Af) J] II MAP-^e-^ + A^) 

j=l,j^i t=±1 j=l 

/ = 2, . . . , n — 1, 
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1 = Y[ e _ 2 (A; n) -A5 n+1) )e_ 2 (Af ) +Af +1) ) 

x J] e^X^ - e^X^ + J] e^A^-Af^e^^+Af) 

i=i j=i,j¥=i 

1 = II e 4 (AS" +1) -Af +1) )e 4 (Al" +1) + Af +1) ) 

J] e_ 2 (Al n+1) - Af >) e_ 2 (A, (n+1) + Aj n) ). (2.44) 
i=i 

In particular the equations for the osp(2\2) open chain are given by 

M< 2 ) 

i=i 

1 = II ^(Af° - e_ 2 (Ap> + Af ) J] ^ - ) e 4 (Af } + A®) (2.45) 

Notice that there is a one-to-one correspondence between the distinguished Dynkin diagrams (see 
fig. 1) and the Bethe Ansatz equations derived for each case. The Bethe Ansatz equations for 
osp(l\2n), osp(2\2n), osp{2m\2), osp{2m + 1|2) can now be compared with the corresponding 
results obtained in [23] for super spin chains with periodic boundaries. Let us point out however 
that we derived explicitly the Bethe Ansatz equations for any osp(M\2n) open spin chain, and we 
expect the corresponding equations for a chain with periodic boundary conditions to be "halved" 
(i.e. half of the factors in the products should be missing) compared to the ones we found. 



2.5 Non-trivial diagonal boundary conditions 

Until now we have derived the Bethe Ansatz equations for trivial boundary conditions, namely 
K~ = K + = 1. We shall now insert non-trivial boundary effects and then rederive the modified 
Bethe Ansatz equations. We choose K~ to be one of the diagonal solutions Dl, D2, D3 found 
in Proposition 3.1 of [3]. We consider, for simplicity but without loss of generality, K + = 1. 
Note that the pseudo- vacuum remains an exact eigenstate after this modification. We rewrite the 
solutions Dl, D2 and D3 of [Sj in a slightly modified notation, which we are going to use from 
now on. 

Dl: The solution Dl can be written in the following form 

K{X) = diag{a, . . . , a, (3, . . . , (5) . (2.46) 

The number of a's is equal to the number of /3's, so that this solution exists only for the osp(2m\2n) 
cases as stated in Proposition 3.1 of ref. |SJ, and 

a(A) = -A + ^, /3(A) = A + z£, (2.47) 
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where £ is the free boundary parameter. 

D2: Solution D2 can be written in the osp(M\2n) case (M > 2) as 



K(X) = diag(l, . . . , 1, a, 1, . . . , 1, P, 1, . . . , 1) (2.48) 



AT -2 



and for the osp(2|2n) case as 



K(\) = diag(a,l,. (2.49) 

2d 



with 



A + z£i A + 

where £i and £ n are the boundary parameters which satisfy the constraint 

6 + en = «-^o- (2.51) 

Obviously, this solution does not exist for the osp(l\2n) superalgebras. 
D3: Solution D3 has the form in the osp(M\2n) case (M ^ 2) 

K(X) = diag(P, . . . , j3, a, . . . , a, j3, . . . , /3, a, . . . , a, {3, . . . , (2.52) 

ji— ni rai+mi M—1m\ m+mi n— ni 

while for the osp(2|2n) case it takes the form 

if (A) = diag(a, a, . . . , a, (3, . . . , (3, a, . . . , a, a). (2.53) 

n\ 2n—2ni m 

The osp(l\2n) case is recovered by taking M = 1 and mi = in (J2.52D . 
Again, a and (3 are given by 

a(A) = -A + z£, /3(A) = A + (2.54) 

where £ = (k + 2m\ — 2ri\ — l)/2 has a fixed value, the integers m\ and n\ being restricted to 
< n\ < n and < mi < m (M = 2m or M = 2m + 1) respectively, ni and mi being neither 
both zero nor taking maximal values simultaneously. 



We now come to the explicit expression of the eigenvalues when K~ is one of the above 
mentioned solutions. We should point out that the dressing functions are related to the bulk 
behaviour of the chain and thus they are form-invariant under changes of boundary conditions. 
Indeed the only modifications in the expression of the eigenvalues (|2.2fij) occur in the gi functions, 
which basically characterise the boundary effects. We call the new gi functions gi. 
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Dl: As already mentioned, the solution Dl can only be applied to osp(2m\2n) with m > 1. 
In this case we have 

&(A) = (-X + i09i(X), l = 0,...,n + m-l 

gi(\) = (\ + i£ + iK,)gi(\), I = n + m, . . . , 2n + 2m - 1 (2.55) 

where (37(A) are given by (|2.24j) - (j2.25j) . The system with such boundaries has a residual symmetry 
sl(m\n), which immediately follows from the structure of the corresponding K matrix. 

D2: We have to separate the cases osp(M\2n) with M^2 and osp(2\2n). In the osp(M\2n) 
case with M^2, one gets 

9i(X) = 9i(X), Z = 0,...,n-1, g n (\) = - 7^7% r^- 9n(X), 

{A + Ki) 

~9i(X) = { ^7YZ 1 r^9iW, l = n + l,...,n + M-2 

[A + t^i) 

(A + i£i - i) (-A + —in + iM — 3i) 

^" +M - l(A) " (A + Z6) (A+^) 9n + M-l(X), 

= (A {X ( ^ n ~^ 9i(X), l = n + M,...,2n + M-l (2.56) 
(A + i£i) (A + <„) 

In the osp(2|2n) case, the formulae are similar to the sp(2n) case: 

= rv i 'n ^2n+l(A) = : — — — — r S-2n+l(A), 

(A + «4i) (A + «4i) (-A - ik + «£i + 1) 

9i(X) = {X ( t^.X %) fli(A), i = l,...,2n (2.57) 
(A + 

The functions (37(A) are given by (|2.24j) - (j2.25j) . The system with such boundaries has a residual 
symmetry osp(M — 2|2n). 

D3: For the D3 solution we find the following modified g functions in the case of osp(M\2n) 
with 

g t (X) = (X + i£)gi(X), l = 0,...,n-m-l 

9i(X) = (—X + i£ — in + ini)gi(X), I — n — tlx, . . . , n + mi — 1 

gi(X) = (X + i---)gi(X), l = n + m 1 ,...,n + m-l 



(\ _|_ {Sl _ L\ 

9m+2u-i-i(X) = .1 I 9i(-X - in), l = 0,...,n + m-l 

{A + l „ + 2 J 



and in the case of osp(2m + l\2n) 



(2.58 



g n +m{X) = (A + z- - -) g n + m (X) (2.59) 
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In the osp{2\2n) case, the formulae become 



gi(X) = (-\ + i£)gi(\), Z = 0,...,m 

gi(X) = (A + i- + -) #(A), / = ni + 1, . . . ,2n- m 

~ 9l (\) = (-A-i/c-ifl ^'I^f ft(A), Z = 2n-m + l,...,2n + l (2.60) 

The system with such boundaries has a residual symmetry osp(M — 2mi\2n — 2n\) @osp(2m\\2n\). 

We now formulate the Bethe Ansatz equations for the general diagonal solutions. The only 
modifications induced on Bethe Ansatz equations are the following for each solution: 

Dl The factor — e 2 ^ 1 +K (A) appears in the LHS of the (n + m) th Bethe equation. 

D2 The factor — e 2 £ 1 _ n (A) appears in the LHS of the n th Bethe equation. 

The factor — e^-n-iW appears in the LHS of the (n + l) th Bethe equation. 

D3 The factor — e 2 £-( n - ni )(A) appears in the LHS of the (n — rii) th Bethe equation. 
The factor — e 25 1 +2ni _ mi _ n (A) appears in the LHS of the (n + mi)" 1 Bethe equation. 

3 Scattering for the osp{\\2n) open spin chain 
3.1 Low lying excitations 

Before we derive explicitly the bulk and boundary scattering amplitudes for the osp(l\2n) case 
we first need to determine the ground state and the low-lying excitations of the model. We recall 
that the energy is derived via the relation H = ^£(A)|a=o- It is given by 

e=~t — - 2 — . (3.i) 

In what follows we write the Bethe Ansatz equations for the ground state and the low-lying 
excitations (holes) of the models under study. Bethe Ansatz equations may in general only be 
solved in the thermodynamic limit iV — > oo. In this limit, a state is described in particular by the 
density functions er'(A) of the parameters A^. 

A. osp(l|2) 

Let us first consider the osp(l|2) case, for which the ground state consists of one filled Dirac sea 
with real strings (all Aj's real). The set of Bethe Ansatz equations for the osp(l|2) case takes the 
form 

M 

eiM^eiMe.iM = l[e 2 (\ i ~ \) e ^ X * + X i) e -i( A * ~ A j) e -i( A * + > ( 3 - 2 ) 

i=i 
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The reason why we study this case separately is basically because we wish to point out the strik- 
ing similarity between the latter Bethe Ansatz equations ()3.2|) and the corresponding equations 
appearing in the study of the 577(3) open spin chain with "soliton non-preserving" boundary 
conditions This is indeed a remarkable connection, which can presumably be extended to 
open spin chains with "soliton non-preserving" boundary conditions, for higher rank algebras. In 
particular, we expect that for any SU(n) (n odd) chain with "soliton non-preserving" boundary 
conditions the resulting Bethe Ansatz equations will have the same form as in the osp(l\n — 1) 
open spin chain with certain diagonal boundaries. We hope to report on this in detail elsewhere. 

We now study the low-lying excitations, which are holes in the filled Dirac sea. In order to 
convert the sums into integrals, after taking the thermodynamic limit (N — > oo), we employ the 
following approximate relation 

jj E ^ = / V(A) -^E /(A?) - — /(0) (3.3) 



where the correction terms take into account the holes located at values Xf and the halved 
contribution at + . For osp(l\2) in particular i/W = v and A^ = Aj. 

We shall denote by f{uS) the Fourier transform of any function /(A). Once we take the logarithm 
and the derivative of (|3.2j) . we derive the densities from the equation 

t{u)a{u) = ai(o;) + ±F{u) (3.4) 

where aAX) = — — — lne/A) and aAu) = e~^. Moreover 
2n dX 

cosh 

4 ,_i 

Finally equation (|3.4|) can be written as 



cosn — / \ 

K{u) = e'^^-^, F(X) = 02(A) - ai (A) + ^ ((a 2 - a x )(A - A») + (a 2 - 0l )(A + A,))- (3.5) 



a(A)= 2 e (A) + 1$(A), (3.6) 

where 

e(w) = ai(u)K(u), ${lu) =K{uj)F(uj), and K(u) = ^(cu) (3.7) 

In particular the energy e can be written in terms of hyperbolic functions as e(u) = ^^at ■ 

B. osp(l|2n) 

We recall that we can only consider the D3 solution in this case. The ground state consists of 
n filled Dirac seas with real strings. With the help of relation ()3.3|) we derive the densities that 
describe the state with i/W holes in the I sea from the equation 

K{u)a{u) = a{u) + ±F(cu) + ±d(u, (3.8) 
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where we have introduced 



a(X) 



I 2ai(A) \ 




a(X) 



( \ 

a\X) 



(3.9) 



,(0 



V / 

F{X) is a n- vector as well with 

F j {\) = a 1 (A)5 jl -a 1 (A) + a 2 (A) + ^(a 2 (A-Af ) ) + a 2 (A + \ (0 ))^ 

i=l 

(ai(A - A?) + ai (A + A?)) (<^ +1 + ^_!), (j = 1, . . . , n - 1) 
ai(A) - ax (A) + a 2 (A) + ^ ((a 2 - a x )(A - A?) + (a 2 - a x )(A + A?))^ 

i=l 

^^(A-Ai^ + a^A + Af))^ 



F n (A) 



(3.10) 



i=l 



and /C is a n x n matrix with entries given by 



K<ij{w) = (1 + a 2 (u;))<% - a 1 (u)(6 itj+1 + i,j = l,...n- 1, 

Knn-iiu) = /C n _i n (a;) = ai(uj), )C nn (uj) = 1 - ai(u;) + a 2 (a;) (3.11) 

Finally, the n- vector G carries all the explicit dependence on the boundary parameter £ of the D3 
solution 



G J (A) = a 2f _^(A) Sjfi 



Hi 



(3.12) 



where h\ — n — n\. Solving equation ()3.8|) we find the densities a 1 which describe a Bethe Ansatz 
state. The solution of ()3.8j) has the following form 



o{X) = 2 e (A) + I$ (A) + i$ 1 (A, 



(3.13) 



where e and $o,i are ^-vectors (columns) with 



?(u) = K il (u)a 1 (u), %{u) = Y,K ij {u)FK"), = ^UijM&fa £) ■ (3.14) 

7?. = /C _1 and e- 7 is the energy of a hole in the j sea; they are written in terms of hyperbolic 
functions as 



sinh ^min(i,j)^j cosh (n + | — max(z, j)j : 



cosh(ri + |)f sinh ^ 



i,j = l,...,n (3.15) 
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3.2 Scattering 

As already mentioned the main aim here is the derivation of the exact bulk and boundary S- 
matrices. We follow the standard formulation developed by Korepin [30 a , and later by Andrei and 
Destri [3T]. One first implement the so-called quantisation condition, 

(e 2iNpl S - 1)^ = (3.17) 

where p is the momentum of the particle (in our case, the hole) with rapidity A^. For the case 
of v (even) holes in I sea we insert the integrated density ()3.13|) into the quantisation condition 
(13.17}) . We use the dispersion relation 

£ '(A) = ^P'(A) (3.18) 

and the sum rule N / dXa(X) E Z + . We end up with the following expression for the boundary 
Jo 

scattering amplitudes: 



a +l a~ l 



exp{27riV^ * d\(a l (X) - 2e'(A)^| (3.19) 



with 

a~ l (X, = fco(A) 0, «+<(A) = fc (A) (3.20) 

where a +l is realised just as the overall factor in front of the unit matrix at the left boundary 
(recall that K + = 1, whereas K~ is given by the solution D3). Moreover, 

fco(Ai) = exp {ztt J 'dA^A)}, h(Xl £)=exp{2inj * d\${(\, £)}, (3.21) 

with <3> z given by (J3.14|) . (|3.1UJI . We finally restrict ourselves to I = 1 in the first sea and we write 
the latter expression in term of the Fourier transform of ([3.14)1 . 

*b(A) = exp { - - J _<|i( w )e~- A }, h(X, = exp { - J — fle"* 1 *}. (3.22) 

Let us discuss first the form of the exact bulk ^-matrix. It is easy to compute the scattering 
amplitude between two holes in the first sea. The bulk scattering amplitude comes from the 
contribution of the terms of $ x given by eqs. (|3.10|) . ()3.14|) . (|3.15|) . with argument A ± Xj. After 
some algebra and using the following identity 

^ 7^37 (3-23) 



2 J u coshf T(^ 
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we conclude that the hole-hole scattering amplitude is given by the expression 

tanvr(l^) r(=%) r(^ + 1) r(=^±i) r(£±i + §) 

°^ ' ~ 7Z^Z7JE+T\ TV -jA \ TV iA i 1 \ TV iA+1 \ T-y -iA+1 i l\^ K 6 -^) 

As a consistency check we compute one further eigenvalue of the S-matrix. In particular if one 
considers the state with two holes in the first Dirac sea, and a two-string located at the midpoint 
of the two holes, the corresponding eigenvalue is given by (see also [32]) 

S b {\) = e 1 (X)S (X). (3.25) 

The explicit bulk S'-matrix, which is a solution of the super Yang-Baxter equation has the following 
structure 

S ^ = r./fun ^ + + i(A + %k)P - iXQ). (3.26) 

We now give the expressions for the boundary S'-matrix, which follow from (|3.21|) , ()3.22|) , and the 
duplication formula for the V function 

2 2x - 1 T(x + -)T(x) =7rir(2x). (3.27) 
The overall factor ko, (|3.21j) . is given by 



r( iX ) r( -»a , 3\ r( iX -i I i- 2 1 ) T( ~ iX -i i i- ±1 

1 \ 2n+l' 1 y2n+l_jJ_ V 2n+1 ^ 2(2n+l) _ 4/ 1 V2n+1 ^ 2(2n+l) ^ 2/ 
-2n+l ^ 1 ^2n+l l) ^(in+T + 2(2n+l) 4) ^(in+T + 2(2n+l) 2' 



0\ / 0^ F('^*A \. tV iA 1 3\. tV — iA 1 1 1 3\ tV iA I 1 \ 

1 V2n 

where 



y /,v _ Sm7r (2^TT + 2(2^11 ~ j) S1I17r (2iTT ~ 2(2n+l) + |j ^^^^TT + 4(2n+l) + ij /o 9q \ 
bLUn \2n+l 2(2n+l) 4/ b111 ^ 2n+l "T 2(2n+l) 2/ LdJ17r v2n+l 4(2n+l) 4/ 

The £ dependent part for the D3 solution fci, (j3.21j) reads 

p( »A I g' I 1\ p/ -iA I \ p/ -iA I g'-ni , In p/ iA , g'-ni , i\ 

1 n £) — V 2n+1 2n+l 2> ^2n+l 2n+lV V2n+1 2n+l 2^ V 2n+l 2n+l ^ /g og\ 

1 ' ~~ TV ~ iA 1 £' 1 In pf iA , pf iA 1 g'-ni 1 1\ p/ -iA , g'-ni \ , inV / 

V2n.+1 ~ 2n+l ' 2> V2n+1 ' 2n+l/ ^2n+l ~ 2n+l ~ 2/ \2n+l ~ 2n+l > ~ > 

where ^' = ^ — | has the fixed value found for D3, so that the boundary S-matrix satisfies the 
reflection equation. Note that our solutions include the necessary CDD factors both for the bulk 
and boundary matrices. 
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